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cn ; Abstract 

o: 

<^ . Polynomials in differentiation operators are considered. The Darboux transformations 

^^ I covariance determines non-Abelian entries to form the coefficients of the polynomials. 

.^ ■ Joint covariance of a pair of such polynomials (Lax pair) as a function of one-field is stud- 

^^ . led. Methodically, the transforms of the coefficients are equalized to Frechet derivatives 

r~| i (first term of the Taylor series on prolonged space) to establish the operator forms. In 

Qh' the commutative (Abelian) case that results in binary Bell (Faa de Bruno) differential 

r^ ■ polynomials having natural bilinear representation. The example of generalized Boussi- 

j^ ! nesq equation is studied, the chain equations for the case are derived. A set of integrable 

^ I non-commutative potentials and hence nonlinear equations is constructed altogether with 

:• ■ explicit dressing formulas. 

^ . 1 Introduction. 

Investigations of general DT theory in the case of differential operators 

L = J2akd^ (1.1) 

with non-commutative coefficients was launched by papers of Matveev pp. The proof of a 
general covariance of the equation 

i^t = Lijj (1.2) 

with respect to the classic Darboux transformation 

^[l\=^'-ai; (1.3) 

incorporates the auxiliary relation 

at = dr + [r,a], , . 



where Bn are differential Bell (Faa de Bruno [5,) polynomials fU\. The relation ()1.4|) generalize 
so-called Miura map and became the identity when a = 4>'(p^^, is a solution of the equation 
(0321). 

Such combinations may be used in the Lax representation constructions for nonlinear prob- 
lems. It opens the way to produce wide classes of solutions of the nonlinear problem. This 
known approach j^ intensely develops nowadays. Examples of non-Abelian and non-local equa- 
tions, integrable by DT was considered in |2j, jlj. Some of them were reviewed and developed 
in the book [2j. The general approach was developed for difference operator (see again the 
papers cited above) and recently generalized for wide class of polynomials of automorphism on 
a differential ring [TUj. 

The choice of the jointly covariant combinations introduces addition problems of the appro- 
priate choice of potentials on which the polynomial coefficients depend [15J. This problem was 
recently discussed in ^^l; where a method of the conditions account was developed. Covari- 
ant combinations of (generalized) derivatives and potentials may be hence classified for linear 
problems. The task is intimately connected with other reduction problems of Lax operator rep- 
resentation of nonlinear equations. In two words, having the general statement about covariant 
form of a linear polynomial differential operator that determine transformation formulas for 
coefficients (Darboux theorem and its Matveev's generalizations), the consistency between two 
such formulas need the special definition of potentials. For example, the second order scalar 
differential operator has the only place for a potential and the covariance generate the classic 
Darboux transformation for it. The place that this potential take in the second equation of the 
KdV Lax pair needs special investigation. 

In scalar case such one-potential construction have been studied in [13] and developed for 
higher KdV and KP equations ^1] . It was found that the result is conveniently written via such 
combinations of differentiation operator and exponential functions of the potential as Binary 
Bell Polynomials (BBP) [llj. The principle and result is reproduced and developed in the Sec. 
1.2 of this paper to give more explanation. 

The whole construction in general (nonAbelian) case could be similar, a bit more compli- 
cated, but much more rich and promising. The theory could contain two ingredients. 

i)The first one would be non-Abelian Hirota construction in the terms of the mentioned 
binary Bell polynomials JT] . On the level of general expressions and applications some obstacle 
appears, e.g. an extension of addition formulas ^T]) to the nonAbelian case. 

ii) The second way relates to some generalized polynomials that could be produced as 
covariant combinations of operators with a faith that observations from Abelian theory could 
be generalized. Namely the case we would discuss in this paper. 

Even for the minimal (first order in the D-operator) examples one arrives to the operator ZS 
problems that contain many interesting integrable models. It is seen already from the point of 
view of symmetry classification [H]. So, the link to DT covariance approach allows to hope for a 
realization of the main purpose - construction of covariant polynomials and their enumeration 
and use in the soliton equations theory. 

We would begin from the example, using notations from quantum mechanics to emphasize 
the non-Abelian nature of the consideration. The operators p and H could play the roles of 
density matrices and Hamiltonians, respectively, but one also can think of them as just some 
operators without any particular quantum mechanical connotations. The approach establishes 
the covariance with respect to bDT [?] of rather general Lax system for the equation 

-ipt = [H,h{p)], 



where h{p) - analytic function, in some sense -"Abelian", the function to be defined by Taylor 
series 1201 • More exactly it is shown that the following statement takes place 

Theorem. Assume {ip\ and \ip) are solutions of the following (direct) equations 

Zu{x\ = {x\{p-^H), 

-i{Xt\ = -{x\Hp), 

V 

and \lp) \lp) is for conjugate 'pair Here p, H are operators left- acting on a "bra" vectors {ip\ 
associated with an element of a Hilbert space and transforms (bDT) complex numbers X, z\ are 
independent oft. and {ipi\, pi, h{p)i are defined by 

{H = m{i + ^p), (1-5) 



p — X 



Pi = TpT-\ h,{p) = Th{p)T-\ T = (l + ^P), (1.6) 

where P = \{p) {x\(p) {x\ ■ Then the pairs are covariant: 

z\{i)i\ = (V'iKpi - ^H), -i{^i\ = -(^i|/ii(p). 

The Lax-pair representation and Darboux covariance properties of this equation have been 
established in [20] • The cases /(p) = ip^ and /(p) = ip~^ were considered in |S|. A step to 
further generalizations for essentially non- Abelian functions, e.g. h{X) = XA+AX, [A, X] ^ 0, 
is studied in ^7] . The case is the development of the matrix representation of the Euler top 
model [211. 

This example of the theory is more close to the spirit of the sec. 2.4 and more achievements 
are demonstrated in [22l, where abundant set of integrable equations is listed. The list is in 
a partial correspondence with [HI, and give the usual for the DT technique link to solutions 
via the iteration procedures or dressing chains. The results show how the "true" non- Abelian 
functions appear in the context of the covariance conditions application. On the way of this 
specification we use the notion that is similar to one for automorphic function. 

2 Covariance principle equations 

2.1 One-field Lax pair for Abelian case. Covariance equations 

First we would reproduce the "Abelian" scheme, generalizing the study of the example of the 
Boussinesq equation ^21- To start with the search we should fix the number of fields. Let 
us consider the third order operator ()1.1|) with coefficients bk,k = 0,1,2,3, reserving a^ for 
the second operator in a Lax pair. Suppose, both operators depend on the only potential 
function w. We would restrict ourselves by the case of bix = and the special choice 63 = 1 , 
61 = b{w), bo = G{w). The problem we consider now may be formulated as follows: To find 
restrictions on the coefficients b2{t),b{w,t),G{w) compatible with DT transformations rules of 
the potential function w induced by DT for bi. The standard DT ()4.45|) for the third order 
operator coefficients (Matveev generalization |1J yields (we denote df = /') 

62[l] = 62 + &3, (2.7) 



bi[l] = h + b'^ + 3ha' , (2.8) 

60 [1] = bo + b[ + ab'^ + 3bs{aa' + a"), (2.9) 

having in mind that the "elder" coefficient 63 does not transform. Note also, that 63 = yields 
invariance of the coefficient 62 • 

The general idea of DT form-invariance may be realized considering the coefficients trans- 
forms to be consistent with respect to the fixed transform of w. Generalizing the analysis of 
the third order operator transformation [T2j, one arrives at the equations for the functions 
62(^)5 b{w, t), G{w). The covariance of the spectral equation 

bsi^^^^ + b2{t)tlJ^x + b{w, t)i)^ + G{w)i) = Xi) (2.10) 

may be considered separately, that leads to the link between 6, only. We study the problem 
of the ()2.10|) in the context of Lax representation for some nonlinear equation, hence the 
covariance of the second Lax equation is taken into account from the very beginning. We name 
such principle as the "principle of joint covariance" fH]- The second (evolution) equation of 
the case is : 

i)t = a2{t)i'xx + ai{t)i^x + wi), (2.11) 

with the operator in the r.h.s. having again the form of ()1.1|) . We do not consider here a 
dependence of aj,6j on x for the sake of brevity, leaving this interesting question to the next 
paper. 

If one consider the L and A operators of the form (1), specified in equations ()2.1()|1 and 
()2.1H1 . as the Lax pair equations, the DT of w implied by the covariance of ()2.1H) . should be 
compatible with DT formulas of both coefficients of ()2.10|) depending on the only variable w. 

a2[l] = a2 = a{x,t), 

ai [1] = ai (x, t) + Da{x, t) 

ao[l] =w[l]=w + a\ + 2a2(r' + aa'^ (2.12) 

Next important relations being in fact the identities in the DT transformation theory |6j, 
see the introduction, are the particular cases of the generalized Miura map. p.4|l : 

at = [a2{a'^ + aa:) + aia + w]x (2-13) 

for the problem flT^ and, for the (jTTTl|l 

cr^ + SaxCr + (Jxx + b{w, t)a + G{w) = const; (2-14) 

is a solution of both Lax equations. Suppose now that the coefficients of the operators are 
analytical functions of w together with its derivatives (or integrals) with respect to x (such 
functions are named functions on prolonged space [IE])- For the coefficient G it means 

G = G{d'^w, w, Wx, ...d'^wt, wt, Wt^, ...), (2.15) 



The covariance condition is obtained for the Frechet derivative (FD) of the function G on the 
prolonged space, or the first terms of Taylor series for ()2.15|) , read 

G{w + a[ + 2a2a' + aa'^) = G{w)^G^^{a[ + 2a2a' ^-aa'^)' + G9~i^Xau + 2a2d'^{a[)+d'^{a) + ..., 

(2.16) 
where we show only terms of further importance; the expression simplifies if 02 does not depend 
on X. Quite similar condition one have for thefei = b{w,t), with which we would start. In the 
analogy with the expressions ()2.16l2.8p one obtains 

b'2 + 363O"' = bw{a[ + 2a2<j' + cra'2) + &«,'(% + 2a20"' + cra'g)'.... (2-17) 

This equation we name the (first) "joint covariance equation" that guarantee the consistency 
between transformations of the coefficients of the Lax pair ()2.11|1 . ()2.10|1 . In the frame of our 
choice 02 = 0, the equation simplifies 

'^2 — ^w(^li 

or 

b^ = 3b3/2a2 , ^ 

b'^ = a[3b3/2a2. ^^^ 

So, if one wants to save the form of the standard DT for the variable w (potential) the sim- 
ple comparison of both transformation formulas gives for b(w) the following connection (with 
arbitrary function) 

b{w,t) = 3w/2 + a{t). (2.20) 

Equalizing the expansion ()2.16|) with the transform of the 60 = G{w) yields: 

b[ + ab'2 + 863(^72 + <y'y = G^M + 2a2a' + a^)' + G^-i^Jait + 2d-\a2a[) + d-\aa'2)t] + ... 

(2.21) 
This second "joint covariance equation" also simplifies when a'g = 0: 

863^7202 + 863(5" Vi - wy/a2 + 863^72 = G^.(2a2(T')' + Ga~i^,[2a2at] + ... (2.22) 

when (I2.2()|l is accounted. Note, that the "Miura" 1)2.1111 linearizes the FD with respect to a. 
Finally, 

Gw:, = 863/202 ^2 23") 

Gd-iwt = 863/202 

Compare with the formula w[l] = w + 2a x- The transformation for the potential w follows 
from the last equation of the system ()2.9|) . i.e., 

G\\\ = G + 3wxl2 + 8(^2/2 + a,).. (2.24) 

Such equation determine the functional dependence of G(u) and we would name such equations 
as joint covariance equations. We see that further analysis is necessary due to such constraint 
(reduction) existence. Then, similar to the case of KdV equation, the covariant equations 
(12.1012.111) 

O'lt + 0-20-1 + aia[ = 0, (2.25) 



which get the form of the Burgers equation after ()2.19|) account. Finally the "lower" coefficient 
of the third order operator is expressed by 

G{w, t) = 3b3W^/2a2 + 3b3a[d'^w/2{a2Y + 3b3d-^Wt/2al (2.26) 

Statement 1 The expressions ()2.1H I2.10| 12.201 I2.26|l define the covariant Lax pair when the 
constraints ()2.251 I2.19|l are valid. 

Such equation for ()2.1H1 compare with Riccati equation (stationary version) for the second 
order spectral problem corresponding to KdV provided by the name of Miura equation. If one 
would use the equation ()3.30|) in ()2.24|) . the time-derivative of w appear. Moreover, the further 
analysis shows that the case we study need to widen the functional dependence in u, namely 
we should include not only derivatives of w with respect to x, but integrals (inverse derivatives) 
as arguments of the potential. Let us introduce analytical function G denoting the coefficient 

bo. 

The DT transform of G after substitution of ()2.24|) gives 

G{w) + 3wj2 + 3{-at - w^)/2 + 3a^j2. (2.27) 

Equahzing (IT77I) and ^HM yields 

Gw^ = 3/4; Ga-i^i = —3/4. 

That leads to the exact form of the Lax pair for the Boussinesq equation from [7j for the choice 
of a = -3/4. 

Remark 1. We cut the Frechet differential formulas on the level that is necessary for the 
minimal flows. The account of higher terms leads to the whole hierarchy [?]. 

Remark 2. We cut the Frechet differential formulas on the level that is necessary for the 
minimal flows. The account of higher terms lead to higher flows (higher KdV, for example) 



3 The solitons of the generalized Boussinesq and dress- 
ing chain for the Boussinesq equation 

To produce a simplest soliton solution of the generalized Boussinesq equation, it is enough to 
start from zero potential in the Lax pair equations (j2. 1012.111) . 

bsi^xxx + b2(t)tp^^ + a{t)ij^ = X^j .g2gN 

Ipt = a2{t)lp:rx + a'l{t)lpx 

The seed solution should satisfy both Lax equations (|3.28|) The dressing formula for the zero 
seed potential ()4.45|) is standard and includes this only function. 



w. 



a[ + 2a2<j' + cra'2 = a[ + 2a2 log^^ 0(x, t) (3.29) 



Going to the dressing chain, we use the method from [Tz]. We would restrict ourselves to 
the case of constant 62 = 0, 63 = 1 , 61 = 6, 60 = u. The general construction is quite similar, 
the covariant equations ()2.11I2.10|) are accompanied by the following equation 

(rt = -{a'^ + a^)^-w^ (3.30) 



for the problem ()2.1H) and 

a^ + ?)(Tx(y + (Txx + ha + G = const, (3.31) 

for ()2.10j) . see (J2.13J) . compare with ()2.14|1 that was Riccati equation (stationary version) for 
the second order spectral problem corresponding to KdV. If one would use the equation ()3.3()j) 
in ()2.24|1 . the time-derivative of w appear. 

Namely the "Miura" equations ()3. 301 IH3T] together with the DT formula 11.51 

Wn+i = Wn + lln^an (3.32) 

form the basis to produce the DT dressing chain equations. 
We express the iterated potential Wn from 13.301 

- (^nt - (0"„ + Crnx)x = Wnx (3.33) 

and substitute it into 13.321 The first dressing chain equation 

(Jn+i,t - (^nt = {(^l+i + CT^+i)' - {crl - (y'J- (3.34) 

Next chain equation is obtained when one plugs the potential from ()3.33|) to the iterated (??) 

o\ + 3o-^o-„ + cr^ + (-3u„/2 + a)a„ + -3m'„/4 + W^Unt = c„. (3.35) 

4 Non-Abelian case. Zakharov-Shabat (ZS) problem. 

4.1 Compatibility condition. 

In the case a'g = by which we have restricted ourselves, the Lax system ()2.11I2.10|) produces 
the following compatibility conditions: 

2036^ = 363a'2, 

bst = ci2&3 + 202^2 + ciib'3 — 36302 — 3b3a[ — 26202 

b2t = 0262 + 2a26'i + ai6'2 - b^a'^ - 6202 - ^i«2 - S^sa'/ - 2b2a[ - ^b-^a'^ (4.36) 

bit = '^261 + aib[ — b^a'C — 62^1 — bia[ — Sb^a'^ — 262^0 + 20269 

bot = cti^o - 61^0 + ^260 - ^aflo - ^sOq" 

In the particular case of 03 = we have at once the direct corollary of the first of equalities 
()4.36|) 63 = 0; in the rest of the equations the restriction ()4.36p is taken into account. 

4.2 Compatibility conditions for two general ZS problems 

Let us list equations that appear as compatibility condition of two first order operators with 
coefficients to form a nonabelian set. We recall these equation just to examine the complete 
set of them from the DT invariance point of view. The pair we study is 

iJt = (flo + aiD)ij, (4.37) 

^|Jy = (60 + biD)ij, (4.38) 



where ai,6j are functions of x,y,t and let D be a differentiation with respect to x operator. 
The equations of compatibihty are 

[6i,ai = 0] (4.39) 

aiy - bit + [ao, h] + [ai, bo] + aib[ - bia[ = 0. (4.40) 

aoy - bot + [ao, bo] + aib^ - bia^ = 0. (4.41) 

The " '" denotes the derivative with respect to x. 

If the coefficients ai, bi do not depend on x, y, t, the first two equations simphfy. 

[&i,ai = 0] (4.42) 

[ao,6i] + [ai,6o] = 0. (4.43) 

Analyzing the conditions, one sees, that the direct proof of the heredity of ()4.40|) is a corollary 
of Jacoby identity. 

If one changes the differentiation operator D in (??) to the shift operator T, the equations 
are changing as follows 

In the next generalization jTUj, the operator T can be considered as automorphism in a 
differential ring. Some direct generalization of (??) is achieved, if the operator D is considered 
as abstract differentiation. See. for example, |T^ where Z) is a commutator. 

4.3 Covariance equations 

Let us change notations 

^t = {J + uD)^, (4.44) 

where the operator J does not depend on x, y, t and the potential ao = u = u{x, y, t) is a function 
of all variables. The operator D now is a differentiation by x. The transformed potential 

u = u+[J,a], (4.45) 

the a = 0x0~^ is defined by the same formula as before. The covariance of the operator in 
()4.37|) follows from general transformations of the coefficients of a polynomial [7j. The operator 
J does not transform. 

Suppose the second operator of a Lax pair has the same form, but with different entries and 
derivatives. 

iljy = {Y + wD)iIj, (4.46) 

where the potential w = F[u) is a function of the potential of the first ()4.44|) equation. The 
principle of joint covariance ^S] hence reads: 

w = w+[Y,a] = F{u+[J,<T]), (4.47) 

with the direct corollary 

F{u) + [r, a] = F{u + [J, a]). (4.48) 

8 



It implies the same functional dependence of the coefficient F on m before and after transfor- 
mation. So, the equation ()4.48|) defines the function F(u), we shall name this equation as joint 
covariance equation. In the case of abelian algebra we used the Taylor series (generalized by 
use of a Frechet derivative) to determine the function. Now some generalization is necessary. 
Let us make some general remarks. 

A class operator- valued function F(u) of an operator u in Banach space may be considered 
as a generalized Taylor series with coefficients that are expressed in terms of Frechet derivatives 
[?]. In a sense of the space norm the linear in u part of the series approximate the function 

F{u) = F{0) + F'{0)u + ... 

The representation is not unique and the similar expression 

F{u) =F{0) + uF'{0) + ... 

may be introduced (some fundamentals about the definitions are given in Appendix). Both 
expressions however are not Hermitian, hence not siutable of a majority of physical models. It 
means, that the class or such operators is too restrictive. To explain what we have in mind, let 
us consider examples. 

4.4 Important example 

From a point of view of the physical modelling the following approximation 

F{u) = H+u + uH 

is preferable for the important class of Hermitian operators. The case of the Hermitian H is in- 
cluded. Such models could be applied to quantum theories: introduction of this approximation 
is similar to "phi in quadro" (Landau-Ginzburg) model []. So, let it be 

F{u) =Hu + uH, (4.49) 

by a direct calculation in ()4.48|) one arrives at the equality 

[Y,a]=H[J,a] + [J,(T]H. (4.50) 

The obvious choice for arbitrary a is Y = H^, J = H. 

The compatibility conditions (?? for ()4.44|) and ()4.46p . when a is Y = H^, J = H yields 

[J, J'] = 
[J,u]H + {u[J,H]}+[u,H]J=0 (4.51) 

Uy — Hut — UtH + [u, H]u + u[u, H] + JHux + JuxH + HJux = 



So this important case of [H,J]=0 gives 



[J,u]H+[u,H]J = 
Uy — Hut — UtH + [m, H]u + u[u, H] + JHux + Ju^H + HJu^ = 



(4.52) 



Both condition are covariant, the first of them is a DT- hereditary constraint, that may be 
checked directly. We hence obtain a class of DT-integrable equations/Lax pairs. 

Uy - {H.u}t + [^i^ H] + JHux + Ju^H + HJu^ = 0, (4.53) 



It is seen that in the case J=H, the first of equations ()4.52j) is vahd identically, 
if the potential does not depend on t it is reduced to the next equation: 

Uy + [u^, J] + J'^u^ + Ju^J + J^u^ = 0, (4.54) 

and x-independence yields the generahzed Euler top equations 

Uy + [m^ J] = 0, (4.55) 

which Lax pair ()4.44l4.46p with Y = J^ was found by Manakov fIT]. So, we used the compati- 
bility condition to find the form of integrable equation and reduction tracing the simplifications 
appearing for the subclasses of covariant potentials. While doing this we also check the invari- 
ance of the equation and heredity of te constraints. 

4.5 Covariant combinations of symmetric polynomials 

The next natural example appears if one examine the link ()4.50p . 

P2{H, u) = H\ + HuH + uH^ 

The direct substitution in the covariance and compatibility equations leads to covariant con- 
straint that goes to identity if, Y = H^, J = H. 

It is easy to check more general Y = J^, J = H connection possibility that leads to the 
covariance of 

n 

Pn{H,u) = Y.H^-^uHK 

Such observation was exhibited in [?] . On the way of a further generalization let us consider 

f{H,u)=Hu + uH + S\ + SuS + uS'^ (4.56) 

Plugging ()4.56|) as F{u) = f{H,u) into ()4.48|) . representing Y = AB + CDE yields 

A[B, a] + [A, a]B+CD[E, a]+C[D, a]E+[C, a]DE = H[J, a] + [J, a]H+S'^[J, a]+S[J, a]S+[J, a]S\ 

The last expression turns to identity iiA = B = J = H,C = aH, D = aH, D = aH, S = (3H 
and [a, H] = 0, [/5, H] = with the link a^ = (3'^. 

Statement Darhoux covariance define a class of homogeneous polynomials Pn{H,u), sym- 
metric with respect to cyclic permutations. A linear combination of such polynomials J2n=i PnPniH, u) 
with the coefficients commuting with u, H is also covariant, if the element Y = J2n=i dnH^n + l) 
and ai = Pi = l,a^"''^ = /5""^^,n 7^ 1. A proof could be made by induction that is based 
on homogeneity of the P„ and linearity of the constraints with respect to u. The functions 
Fh{u) = J^o" (^nPHn{u) Satisfy the constraints if the series converges. 

5 Conclusion 

The equation () generalize Boussinesq equation and has direct physical sense for both case of 
waves launched by initial or boundary problem. In the first case t-dependence of coefficients 
means the external conditions varied with time. In the second case the T and x coordinates 

10 



interchange: the coefficients dependence may be interpreted as conditions that are changed 
when a wave propagates (e.g. a bottom slope for the surface waves). The main result of this 
paper is the covariant equation (j4.48|) or its version (J4.5(J|) . A class of potentials from |22] 
contain polynomials PHn{u) and give alternative expressions for it. The linear combinations, 
introduced here could better reproduce physical situation of interest. 
The work is also supported by KBN grant 5P03B 040 20. 
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6 Appendix. Right and left Frechet derivatives 

The notion of a derivative an operator by other one is defined in a Banach space B [?]. Two 
specific features in the case of a operator-function F{u) G Bu E B should be taken into 
account: a norm choice when a limiting procedure is made and the nonabelian character of 
expressions while the differential and difference introduced. 

Definition. Let a Banach space B have a structure of a differential ring. Let F be the 
operator from B to B' defined on the open set of B. The operator is named the left-differentiable 
m. uq E B if there exist a linear restricted operator L{uq), acting also from B to B' with the 
property 

L(mo + h) — L{uq) = L{uo)h + a{uQ, h),\\h\\ ^ 0, (6.57) 

where ||a(Mo, Z?-)!!/!!/?-!! -^ 0. The operator L{uq) = F'{uq) is referred as the operator of the 
(strong) left derivative of the function F{u). The right derivative F'{uo) could be defined by 
the similar expression and conditions, if one changes Lh — > hL in the (??). 
The Gatoux derivative 
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